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The scattering of two-dimensional (2D) massless Dirac electrons is investigated in the presence of 
a random array of circular mass barriers. The inverse momentum relaxation time and the Hall factor 
are calculated and used to obtain parallel and perpendicular resistivity components within linear 
transport theory. We found a non zero perpendicular resistivity component which has opposite 
sign for electrons in the different K and K' valleys. This property can be used for valley filter 
purposes. The total cross-section for scattering on penetrable barriers exhibit resonances due to the 
presence of quasi-bound states in the barriers that show up as sharp gaps in the cross-section while 
for Schrodinger electrons they appear as peaks. 
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I. INTRODUCTION 

Nanostructures have become the system of choice for 
studying transport over the past few yearsi. Start- 
ing with 2D (two dimensional) electron systems at 
the interface of two different materials several decades 
ago^, recently it has shifted to 2D relativistic materials 
e. g. graphene^"— and topological insulators^. In pristine 
graphene the conduction and valence bands touch each 
other in six points of the Brillouin zone and are defined 
by two independent sets of cones commonly called K and 
K' . Near these points the electronic dispersion is linear 
which corresponds to the dispersion of massless relativis- 
tic particles, described by the Dirac- Weyl equation^. 

During the last decades there has been a lot of theo- 
retical and experimental attempts to use the spin of the 
electron as a carrier of informationiS. Graphene in ad- 
dition to the spin of the electron has two more degrees 
of freedom, sublattice pseudospin and valley isospin or 
valley index. In order to scatter an electron from the K 
valley to the K' valley a large transfer of momentum is 
needed. Typical disorder and Coulomb-type of scatter- 
ing is unable to provide this momentum and in such a 
case the valley isospin is a conserved quantum number in 
electronic transport. This allows us to use valley isospin 
as a carrier of information. It was shown that graphene 
nanoribbons with zigzag edges^i can be used as a valley 
filter. Another promising possibility to control the val- 
ley index of electrons is by using line defectsi^. These 
can be formed in graphene when grown on a Nickel sub- 
strate or by using so called mass barriers that can be cre- 
ated by proper arrangement of dopants in the graphene 
sheets. 

The purpose of this work is to apply the above men- 
tioned ideas to electron transport in the presence of circu- 
lar mass barriers. We solve the Dirac electron scattering 
problem on a sharp circular mass barrier and calculated 
the cross-section, the inverse momentum relaxation time 



and the probability for the electron to be reflected in the 
perpendicular direction. In spite of the circular symme- 
try of the scatterers we obtain a non zero perpendicular 
component of resistivity that allows us to separate elec- 
trons with different chirality, or belonging to K and K' 
valleys. The scattering of Dirac electrons by a penetra- 
ble circular mass barrier is influenced by the presence 
of quasi-bound states that results in resonant behavior. 
The obtained results are compared with those for stan- 
dard Schrodinger electrons. 

The paper is organized as follows. In section II we in- 
troduce the problem and illustrate our formalism by con- 
sidering first the more simple problem of Dirac electrons 
scattered by a circular hole. In section III the boundary 
condition is obtained for a circular barrier i. e. a mass 
dot and a mass anti-dot. The bound states in the mass 
dot are calculated in section IV. In section V the scatter- 
ing by an impenetrable circle is considered and the total 
cross section, relaxation time and Hall component of the 
resistivity is obtained. In section VI we repeat this cal- 
culation for a penetrable circle and compared the results 
with the results for standard electrons that are presented 
in Appendix A. Our conclusions are presented in section 
VII. 



II. PROBLEM 

We consider a Dirac electron interacting with circular 
barrier structures shown in Fig. [T] by the shadowed re- 
gions. In the long wave approximation it is described by 
the stationary equation 



{H -Ej-^ = 
with the following Dirac Hamiltonian: 
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Figure 1: (Color online) Layout of the mass barrier profile: 
(a) - dot, (b) - anti-dot, and (c) - penetrable circle, i. e a 
ring shaped barrier. 



This 4x4 matrix Hamiltonian describes low energy exci- 
tations in K and K' valleys. Due to its diagonal form it 
is possible to separate the scattering problems in each K 
or K' valley. We choose the following 2x2 Hamiltonians: 



H 
H' 



<T p + K(r)(T^, 

(T* p K'{r)a^ 



(3a) 
(3b) 



where cr = {ax,(Jy} and az stand for the Pauli matri- 
ces, and K, (r), K'(r) characterizes the mass barrier for 
K and K' electrons, respectively. We use dimensionless 
variables where velocities are measured in Fermi veloc- 
ity unit V p , all coordinates are measured in the radius R 
of the circular scattering barrier, shown in Fig. [1] by the 
solid red line, and the electron energy — in hvp/R units. 
From now on all equations are for K valley electrons, 
except if otherwise specified. 

According to standard scattering theory we present the 
wave function as 
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consisting on the incident wave in the x-direction and 
the scattered part. The scattered part is characterized 
by the two component scattering amplitude 



(5) 



As the incoming wave function part is normalized to unit 
flow density the differential cross-section is equal to the 
radial flow of electrons corresponding to the scattered 
wave function part, namely. 



and two more averages: the inverse electron momentum 
relaxation time (the quantity proportional to the dissi- 
pative component of resistivity) 



7 = / (i(y5(T((^) (1 — cos (ys) 
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and the quantity 



rj = d(p<7{ip) sin i^, 

^0 
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that corresponds to the perpendicular component of re- 
sistivity (or the analog of the Hall component in the case 
with magnetic field). 

Due to the azimuthal symmetry of our problems it is 
convenient to use polar coordinates. Thus, assuming the 
total wave function as 



*(r) 



t/(r) 



we shall solve the set of equations 
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V = Q{E-k)U, (11a) 
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for the wave function components. 

Now we expand the wave function into partial waves 



replacing Eqs. ([TT|) by the following set of two radial 
equations: 



d m + 1 
dr r 



Vm ^ {E - K)Ur. 



(13a) 



d TTl 

) M„i = -{hi + E)v„i. (13b) 

dr r 



These equations are our main instrument when consider- 
ing the problems shown in Fig. [TJ 



III. BOUNDARY CONDITIONS 



a(^) = /+(^) (<T • n) /(^) = e'^AB* + e-''^A*B, (6) 

where n ~ (cos ip, sin ip) is the unit vector in the consid- 
ered direction. 

Besides the above differential cross-section we shall 
consider the total cross-section 



a = / dipa{ip), 
lo 



(7) 



For the sake of simplicity we restrict ourselves to model 
problems with large mass potentials, and replace the po- 
tentials by proper boundary conditions on the electron 
wave functions. In fact, this is the standard way of de- 
scribing low energy scattering. For this purpose we have 
to solve the appropriate equations in the barrier region, 
to apply the standard boundary conditions for both wave 
function components on the solid circles shown in Fig. [TJ 
and to calculate the limit k — oo. We start with the 
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system shown in Fig. [IJc). Thus, we have to solve the 
fohowing approximate set of radial equations: 



dr 
d 

dr 



(14a) 
(14b) 



in the thin shadowed region // which is delimited by two 
circles of radius 1 ± (S ((5 ^ 1). Its solution reads 



(15a) 
(15b) 



Now satisfying the boundary conditions on the two circles 
demarcating the region // we obtain the following set of 
four algebraic equations 



ui{l-5) = Fe-«* + Gc''^ 
vi{l-5) = -Fe-^^^Ge^^ 



uiii{l + 5) = Fe« 



Ge 



— K.(5 



viii{l + 5) = -Fe^-^ + Ge 



— K.(5 



(16a) 
(16b) 
(16c) 
(16d) 



Eliminating the coefficients F and G we obtain in the 
limit 5 <C 1 the boundary conditions 



viii{l) ~ vi{l) = -{uiii{l)+ui{l)}t&n\i{K5), (17a) 
uui{l)-ui{l) = -{viii{l)+vi{l)]t&Txh{K5), (17b) 

connecting wave function components in the regions / 
and ///. 

For a very thin and very high mass barrier (the analog 
of 5-type barrier for Schrodinger electrons, see Sec. \^ 
we take the following limits: 

5 — )• 0, K — ?► oo, tanh(K(5) = P = const, (18) 

what enables us to rewrite Eqs. (fT7|) as 

vui{l)-vi{l) = -P{uni{l) + ui{l)), (19a) 
uui{l)-ui{l) = -P{viu{l) + vi{l)}. (19b) 

These boundary conditions can be formally replaced by 
inserting Dirac (5-functions into Eqs. namely, replac- 
ing those equations by 



771+1 



dr 



= [E-2P5{r-l)]'. 



(20a) 
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if we assume the following rule for calculating the inte- 
grals when the integrand is a product of Dirac (S-function 
and some function f{x) with the discontinuity^^: 

lim /" dx5{x)f{x) = \ {/(+0) + /(-O)} . (21) 



The parameter P characterizes the effective strength of 
the (5- type barrier, and never exceeds the value P = 1 
in contrast to the case of Schrodinger electrons where it 
can take any value. A second major difference is that the 
wave function components u{r) and v{r) are discontinu- 
ous at the position of the (5-function (but the probabil- 
ity density is continuous) while for Schrodinger electrons 
the wave function is continuous but the derivative of the 
wave function is discontinuous in that position. This is 
a consequense of the fact that the Dirac- Weyl equation 
is a first order diferential equation while the Schrodinger 
equation is second order. 

The obtained boundary conditions for the general case 
shown in Fig. [ijc) enables us to construct analogous 
boundary conditions for the two other cases shown by 
Figs. [IJa,b). So, in the case of Fig.[IJa) we assume 

K = cx), P=l, u///(l) =z;/7/(l) = 0, (22) 

and rewrite the boundary conditions given by Eqs. (jl9p 

as 



Ui{l)=Vl{l). 



(23) 



This is the boundary condition for the quantum dot, sur- 
rounded by an infinite mass barrier. 

In the case of Fig. [TJb) an analogous reasoning leads 
to the following boundary condition: 



uj{l) = -vi{l), 



(24) 



that we shall use for describing Dirac electron scattering 
by a hard wall anti-dot. 

The obtained boundary conditions enables us to ne- 
glect the K, terms in Eqs. ([T^ and solve the Dirac equa- 
tions for free electrons 



771 + 1 



dr 



Vm = EUr. 



dm 

I u„i = -Fv„ 

dr r 



(25a) 
(25b) 



in regions / and /// separately. 



IV. BOUND STATES IN THE DOT 

The most simple problem is the one of a quantum dot 
shown in Fig. [IJa). In this case the solution of Eqs. (pS)) 
in region / has to be finite at r = and is given by Bessel 
functions 



u{r) = FJrniEr), 
w(r) = FJ„,+i{Er). 



(26a) 
(26b) 

Satisfying boundary condition (j23p we immediately ar- 
rive at the algebraic equation 

JraiE) - Jr,-,+l{E) ^ (27) 

that determines the energy of the bound states. 
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V. SCATTERING BY AN IMPENETRABLE 
CIRCLE 

Now we shall consider our main problems, namely, the 
scattering of Dirac electrons by circular mass barriers. 
We shall start with the case presented in Fig. [TJb) — the 
scattering by an impenetrable circle (or by a hard wall 
anti-dot). According to standard scattering theory we 
have to construct the total wave function (jTU]) solving the 
Dirac equations for free electrons ((25|) in the outer region 
/ that satisfies the boundary condition ((24|). to exclude 
the incoming part of the wave function, and calculate the 
cross-section by means of Eq. 

So, the solution of Eqs. ([25]) in the outer region / reads 



u„j(r) = w,n [J,niEr) cos S„i + Y„i{Er) sin S„i] , (28a) 
v,n(r) = [J,n+i{Er) COS 6,n + Y,n+i{Er) sin (5„J . (28b) 

where the symbols J,„ and Y„i stand for Bessel and Neu- 
mann functions, respectively. 

Now satisfying the boundary conditions ((24)) for any 
partial wave function harmonic individually we obtain 
the expansion coefficients sin Sm and cos 5„i through the 
so called phase shifts 



tan Sri 



Jm{E) + J,n+l{E) 
'Y,n{E)+Y,n+i{E)' 



(29) 



Usually the exclusion of the incoming plane wave from 
the total wave function (|12p is done in the asymptotic 
region where Er ^ 1 . Here w e use th e asymptotic of the 
Bessel functions J„i{Er) w yJ^jnEr cos A,„, Y^iEr) « 
sin A,„ where 



A,„ = Er - -^m- -, 



(30) 



which allows us to write the total wave function compo- 
nents as 



C/(r) 



nkr 



w™e""'^cos(A^-<5„), (31a) 



Now in order to compensate the incoming term in the 
total wave function (fT^ by the second term of the last 
line in Eq. ([32]) we have to take 



1 
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(33) 



It is remarkable that this choice makes the above men- 
tioned compensation in both scattering amplitude com- 
ponents A and B, but not in every pair of radial compo- 
nents Um and Vm separately. This is related to the fact 
that the angular momentum is not conserved quantity 
but that the Dirac- Weyl Hamiltonian commutes with op- 
erator consisting of the angular momentum plus isospin. 
So, inserting the above expression of Wm into Eq. ((3T|) 
and subtracting the incoming wave function part given 
by Eq. (|32p we arrive after laborious but straightforward 
calculations at the following expression for the compo- 
nents of the scattering amplitude ([5]): 

^ = £ e""^e-'^-sin(5™, (34a) 

?n — — oo 
r— oo 

^ ^ p-3«7r/4gjv / e™'^e-'*" sinJ,„. (34b) 
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Inserting them into Eq. ^ we obtain the following dif- 
ferential cross-section: 

= ^ E e*[^"""'^^"^'""'"''^]sin5™sin,5„,. 

?n,m' — — OO 

(35) 

We have to keep in mind that the derivation of this 
cross-section was performed with the Hamiltonian pap , 
which is valid for electrons in the K valley. 

To obtain the results for K' valley electrons we have 
to repeat the procedure starting with Hamiltonian (j3b|) . 
Fortunately, it leads to the same set of differential equa- 
tions (fTTj) with a single change U ^ V. Further, we 
replace Eq. ([T^ by 




oo 

E 




(36) 



i(m-\-l)^p 



In this asymptotic region the incoming plane wave can 
be presented as 



m— — oo 



and arrive at the same radial equation set (|13[) as was 
obtained for the K valley case. Consequently, all equa- 
cos(A„i-|_i — Sm)- (31b) tions, including the boundary condition, remain the same 
for the K' valley as well. The equation for the differential 
cross-section for the K' valley, however, differs. In this 
case the changes m — > —m and m' — > — m' have to be 
performed in the argument of the exponent leaving the 
same indexes of phase shifts Sm and Sm' ■ 

Now inserting the obtained differential cross-section 
([55)) into Eq. ([7]) we obtain the total cross-section 



nEr 



J2 *"e'"'^cosA^ 



(32) 



E 



Sm,m'e^^^'^' ^""^ sin Sm sin 5 m' 



m— — oo 



m,m — — oo 



— y 



(37) 
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Note the above mentioned change m — >■ —m and m' — > 
—m' in the exponent argument now appears as the same 
change in the argument of the Kronecker symbol Sm.vi'- 
It is evident that this change does not influence the total 
cross-section, and consequently, it is the same for both 
K and K' valley electrons. 



Taking into account Eq. the partial cross-section 
contribution to the total cross-section can be presented 
as 



sin^ dr, 



[JmiE) + Jm+l{E)] 



(38) 



[JmiE) + Jm+l{E)f + [YmiE)+Yra+l{E)f' 



what enables us to calculate the scattering cross-section 
directly. 

The energy dependence of the total cross-section where 
the sum is restricted by the value M [\m\ ^ M) is shown 
in Fig. [21 We see a rather good convergence at low en- 
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Figure 2: (Color online) The energy dependence of the partial 
sums (J A/ contributing to the total cross-section. The italic 
numbers on the curves indicate the number M. The black 
curve for M = 15 corresponds to the convergent result. 



ergies where three terms (i. e. 
sufficient. 



0, ±1) are already 



The oscillating behavior of the energy dependence of 
the partial sum follows from the same behavior of the 
separate terms in Eq. p7|) that can be easily explained 
calculating the asymptotic of the phase shifts. Indeed, 
using the asymptotic of the Bessel functions and replac- 



pression for the scattering phase (P^l) : 

J.,n{E) + Jrn+l{E) COS A, 



tan S„ 



COS A,„+i 



Y,n{E) 

cosA„, -f- 



^Y,n+l{E) 

sin A„, 



sin A, 
sin (A,„ 



f sin A 
n/4) 



m+l 



sin A„i — cos A„i 
-tan(A„i + 7r/4), 



cos (A,„ + 7r/4) 



Thus, in the asymptotic region we have 
sin^ 5m = s'm^{E - TTm/2), 



(39) 



(40) 



what explains the waving behavior of the obtained partial 
contributions to the total cross-section. 

By the way, this simple expression for the scattering 
phase enables us to perform the approximate summation 
in Eq. ([57)1 for large energies which results in the limit 
cross-section criim = 4 as can also be seen clearly in Fig. [5] 
This value is twice larger than the classical value ad = 2 
that can be obtained assuming that relativistic electrons 
are moving along trajectories given by non quantum me- 
chanical equations of motion. This discrepancy is caused 
by the diffraction of the electronic waves when they are 
scattered by hard wall type potentials, and it is inherent 
over scattering by small angles. It is remarkable that rel- 
ativistic electrons exhibit the same feature as Schrodinger 
electrons (see the textbooki^). 

In Appendix A we present similar results for 
Schrodinger electrons (see Fig. [9]) . Note that there are 
some differences in their fc-dependence: 1) in the low en- 
ergy limit the cross-section of the Schrodinger electrons 
diverge logarithmically, while for Dirac electrons it be- 
comes zero, and 2) cr(fc) for Schrodinger electrons is an 
uniform decreasing function of k while for Dirac elec- 
trons it exhibits oscillations in the low energy region. 
Both cross-sections approach the high energy limit from 
above. 

Inserting differential cross-section (|35|) into Eq. ([8]) we 
obtain the inverse momentum relaxation time (or the dis- 
sipative resistivity component) 



2 



^ ^ [2 Sm' ,771+1 ^m' ,771— l] 

7ri.,77l' — — oo 

X e*^^"^' sin 6m sin Sm' 



(41) 



2 

E 



M 



sin^ {5m - 5m+l) ■ 



ing A,„ by A„ 



we have the following asymptotic ex- 



By the way the change m — > — m and m! — > — m' in 
the arguments of the Kronecker symbol does not influ- 
ence the value of the above expression. Consequently, 
the above inverse momentum relaxation time expression 
is the same for both K and K' valley electrons. 

Now inserting the phases obtained by solving Eq. 
into Eq. (|¥T|) we get the result that is shown in Fig. |31 
It's behavior is qualitatively similar to the one for the 
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Figure 3: (Color online) The contribution of the partial sums 
7j\/ to the inverse momentum relaxation time. The italic num- 
bers on the curves indicate the number M where M = 15 
corresponds to the convergent result. 

total cross-section. For large energies it approaches the 
limiting value 7 — 8/3, which we obtain by calculat- 
ing the integral ([8|) with the classical differential cross- 
section, confirming the known fact that the integral ([8]) is 
not sensitive to forward scattering. Thus this relaxation 
time isn't affected by the above mentioned discrepancy 
between the quantum and classical result as it was with 
the total cross-section. 

Note that there is an essential difference in the separate 
contributions to the total cross-section a and the inverse 
momentum relaxation time 7. The partial contributions 
to 7 do not exhibit any oscillating behavior that was 
inherent in the case of a. This is expected from Eq. ([M)) 
where the difference of neighboring phases {5m — 5m+i) 
doesn't depend on energy in the asymptotic region. 

And at last inserting the differential cross-section ([55)) 
into Eq. ^ we obtain the perpendicular (or Hall) com- 
ponent of the conductivity 

2i °° 

V = 7; / , [Sm',m+1 — Sm'.m-l] 

m,rn' — — oc 

X e'(*™'-«")sin(5„,sin5„,' (42) 

4 

= — ^ sinJm sin (5^+1 sin ((5^+1 - (5m) . 

m— — 00 

The result is shown in Fig. |4l Naively we would expect 
that 77 = at zero magnetic field. To our surprise we find 
that 7/ < and that it conserves it's sign as a function 
of E. It means that the mass barrier acts similar as a 
magnetic field. In order to obtain the result for K' valley 
electrons we have to change m ^ —m and m! ^ —m! in 
the arguments of Kronecker symbols in the first line of 
Eq. (H^ . It is evident that due to this change the Hall 



Figure 4: The Hall component of resistivity. 

component of the resistivity 77 changes it's sign. Thus, the 
electrons from different valleys are deflected to opposite 
sides of the sample. There is no net charge build up 
across the sample and thus no Hall voltage. But there 
is a separation of different K and K' valley electrons 
across the sample and thus we can use this effect for 
valley filtering purposes. 

VI. SCATTERING ON A PENETRABLE 
CIRCLE 

Now we turn to our last problem — scattering of Dirac 
electrons by a penetrable circle in order to demonstrate 
how possible quasi-bound states reveal themselves in the 
scattering cross-section. For this purpose we have to 
solve the Dirac equation for free electrons ([25]) in both / 
and /// regions and to apply boundary conditions ([T^ . 
These solutions are given by Eq. (|26)) for the inner region 
/, and by Eq. (|28|) for the outer region ///. Moreover, 
the procedure of the exclusion of the incident exponent is 
the same as it was performed in Sec.|Vl Thus we can im- 
mediately write down Eqs. (p4|) for the components of the 
scattering amplitudes, and use the previous expressions 
for the scattering cross-section p5|37|41|42p . 

The single procedure that should be performed is to 
satisfy the boundary conditions ([T9)) and calculate the 
phase shifts 5m- Inserting into Eq. (fT9| the solutions 
and we obtain the set of two equations: 

CJm.+l + SYm+1 — FJm+1 

= -P {CJm + SYm + FJm} , (43a) 

= -P {CJm+1 + SYm+1 + FJm+l} , (43b) 

where for the sake of shortness we denoted 

C^Wm cos 5m, S^WmSm5m, (44) 
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and omitted the arguments E of all Bessel functions. 

Now excluding coefficient F we obtain a single equa- 
tion, ft can be solved for the tangent of the phase shift, 
and using the expression for the Wronskian of the Bessel 
functions we arrive at 



to the equation 



tan Sri 



where the symbol 



P 



[l-Pf 



(45) 



(46) 



characterizes the impenetrability of the circle. The value 
p = oo corresponds to a completely impenetrable circle, 
i. e. the previously considered scattering on an impen- 
etrable anti-dot, while the value p = corresponds to 
the case of complete penetration, or the absence of any 
scatterer. 

The numerical results for the lowest contribution (m = 
0) to the total cross-section are shown in Fig. [S] for dif- 
ferent p values. The vertical lines indicate the ener- 




Figure 5: (Color online) The partial ctq contribution to the 
total cross-section for different values of the penetrable pa- 
rameter p of the circular scatterer. The curve marked by 
p = oo corresponds to the case of scattering by an impenetra- 
ble circle. 

gies of the bound states of the dot obtained by solving 
Eq. ((27|) as described in Sec. IIVI Although now the dot 
is penetrable and it has no bound states, the correspond- 
ing quasi-bound states reveal themselves as narrow gaps 
close to the maxima of the oscillating partial contribu- 
tion. Note that this is a particular feature of Dirac elec- 
trons. While in the case of Schrodinger electrons the 
quasi-bound states appear as peaks in the cross-section 
(see Figs. [TUl and [TT] in Appendix A). 

According to Eq. it seems that there should be 
one more set of gaps in the partial cross-section, related 



J„,iE) + J.,n+liE) = 0. 



(47) 



But in this case, however, after neglecting the last term 
in the denominator of Eq. (^51) it coincides with the phase 
([2^ obtained for the scattering by the impenetrable cir- 
cle, and in such a way it indicates that the above condi- 
tion defines the flat minimum related to the diffraction 
pattern in the cross-section considered in Sec. |Vl 

These gaps in the partial cross-section reveal them- 
selves in the total cross-section as shown in Fig. [51 
Because of the contribution of the other partial waves 
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Figure 6: (Color online) The total cross-section a for a pen- 
etrable circular scatterer for different p values. The curve 
marked by p = oo corresponds to the case of scattering on an 
impenetrable circle. 

the gaps no longer reach zero as e. g. shown in the case 
of the partial cross-section cto. Notice that they become 
more pronounced and narrower when the parameter p 
increases (when the circle becomes less penetrable). 

In Fig. [7] the results for the inverse momentum relax- 
ation time and in Fig. [5] those for the perpendicular com- 
ponent of the resistivity are presented. We see that al- 
though Eqs. (|^T|) and ([^^ are more sophisticated func- 
tions of the phases Sm and even two neighboring phases 
dm and Sm+i are intermixed, nevertheless the resonant 
behavior, i. e. the negative peaks, is still clearly visible in 
the inverse momentum relaxation time. The behavior of 
peaks in the Hall component of the resistivity, however, 
is more complicated, and exhibits both sharp peaks and 
dips where now the sign of rj can change in small regions 
of energy. 



VII. CONCLUSIONS 

We investigated the scattering of Dirac electrons by 
sharp circular mass barriers, where we studied both hard 
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Figure 7: (Color online) The inverse momentum relaxation 
time calculated according Eq. (|41l) for a penetrable circular 
scatterer. 



0.5 




Figure 8: (Color online) The Hall component 77 calculated 
according Eq. (|42|) for a penetrable circular scatterer. 



wall type anti-dot and circular penetrable scatterer. For 
this purpose the proper boundary conditions for Dirac 
equation were derived, and it was illustrated how it is 
possible to use formally (5-type functions describing rela- 
tivistic systems with high and sharp potentials. 

The differential and total cross-section, the inverse mo- 
mentum relaxation time and the perpendicular (Hall) 
component of resistivity were calculated. The obtained 
results were compared with analogous results for scatter- 
ing of Schrodinger electrons by similar scatterers. 

It was shown that the scattering of Dirac electrons even 
by azimuthal symmetric structures depends on the valley 
index: the K valley electrons are preferentially deflected 
to one side (the Hall component of the resistivity isn't 
zero) while the electrons of the other K' valley are de- 
flected to the other side of the sample. This enables one 



to use this property for valley index filtering in transport 
experiments. 

There is an essential difference in the energy depen- 
dence of the cross-section between Dirac and Schrodinger 
electrons. At small energies the cross-section for Dirac 
electrons tend to zero while those for Schrodinger elec- 
trons diverge logarithmically. This feature of Dirac elec- 
tron is caused by the fact that zero energy for Dirac elec- 
tron actually corresponds to the middle of the half-filled 
band, and not to the bottom of it as in the Schrodinger 
electron case. 

We showed that in the case of Dirac electron scatter- 
ing on a penetrable circle the quasi-bound states reveal 
themselves as sharp gaps in the total cross-section, the 
inverse momentum relaxation time and the Hall compo- 
nent of the resistivity as well. In the case of Schrodinger 
electrons those resonances show up as peaks in the cross- 
sections, and thus their appearance is qualitatively very 
different. 
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Appendix A: Scattering of a Schrodinger electron by 
circular barriers 

In this Appendix we study the scattering of 
Schrodinger electrons by sharp circular potentials shown 
in Fig. [TJ This allows us to compare them with results 
for Dirac electrons considered in the present paper. 

Here scattering is now described by the single com- 
ponent wave function ^'(r) satisfying the stationary 
Schrodinger equation 

(V^ + /c2)*(r) = 0, (Al) 

where the symbol k stands for the electron momentum 
related to it's energy as E = Now we use dimen- 

sionless units defined as follows. The coordinates will be 
measured in the radius of the circular scatterer R, energy 
— in %^ /mB? units, and the electron momentum — in 
h/R units. 

In polar coordinates the wave function ^'(r) is usually 
expanded into a Fourier series like U (r) component of 
Dirac function in Eq. (jl2p with radial components V'm {t) 
satisfying the Bessel equation. That is why all mathe- 
matics is practically the same as used in previous sections 
with a single replacement of the Dirac electron energy E 
by the momentum k of the Schrodinger electron. Con- 
sequently, we obtain the same differential cross-section 
given by Eq. ([SS]), and Eqs. ([371) and (gl]) for the total 
cross-section and inverse momentum relaxation time. 

Nevertheless, there is an essential difference between 
Dirac and Schrodinger electrons which is the different 
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boundary conditions that will lead to different scattering 
phases. 

Thus, in the case of scattering on a circular hard wall 
potential (see Fig. [TJa)) every radial component of the 
wave function has to satisfy the zero boundary condition 

^-^(1) = 0, (A2) 
what leads to the following equation for the phase 

Jjyi{kCl) 



tan 6„ 



r„,.(fca) '' 



(A3) 



instead of Eq. for the Dirac electron. The partial 
contributions and the total cross-section calculated by 
using the above phase equation are shown in Fig. 1^1 
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Figure 9: (Color online) The energy dependence of the 
partial sums aM contributing to the total cross-section for 
Schrodinger electrons. The italic numbers on the curves indi- 
cate the number M. 

In the case of an extremely narrow penetrable circle, 
i. e. a 5-type potential, Eq. (|Aip has to be replaced by 



{V^ + fc2 - pS{r - 1)} ^{r) = 



(A4) 



what leads to the following boundary conditions for the 
radial wave function components on the circle: 

V',„(l + 0) = ^,„(l-0), (A5a) 
C(l + 0) - V4(l - 0) = (A5b) 

This leads to the following scattering phase equation: 



tan(5„ 



y„(fc)-2/p7rJ^(fc)- 



(A6) 



From it we obtain the following partial contribution to 
the total cross-section: 



■ sin^ Sr, 



k{Jlik) + [Y^ik) - 2/p7rJmikWy 



(A7) 



Figure 10: (Color online) The m — contribution to the 
total cross-section for Schrodinger electrons scattered on pen- 
etrable circular potentials, shown by red soUd curve. Green 
dotted curve is the envelop function 4/fc, and the blue dashed 
curve is the result for scattering on impenetrable scatterers. 



The typical contribution(when m = 0) is shown in Fig. 1101 
by red solid curve. Narrow peaks appear close to the po- 
sitions of the bound states of a dot that are defined by the 
equation Jm{k) = 0. This can be also seen from Eq. (jA7P 
which formally is similar to a Lorentzian curve. The top 
of the peak is achieved when the second term in the de- 
nominator (the analog of detuning) in Eq. (jA7|) is zero. 
In the case of small penetrability of the scatterer this can 
be realized if the large parameter p 3> 1 is compensated 
by a small Jm{k) <^ 1 value. But then the contribution 
becomes equal to 4/fc what indicates that the maximum 
of all peaks reach the above envelope function shown by 
green dotted curve in Fig. 1101 One more property of the 
partial contribution follows from the fact that it is rather 
close to the same contribution in the case of scattering by 
impenetrable scatterer shown by the blue dashed curve 
(Eq. (jA6|l converts itself into Eq. (|A3P when p = oo) 
what indicates that the peaks appear at the minima of 
that dashed curve. 



These resonances show up also in the total cross- 
section as seen in Fig. [TTJ The comparison with the blue 
dashed curve calculated for p = oo indicates clearly that 
the quasi-bound states appear in the total cross-section 
as positive peaks. 
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Figure 11: (Color online) The same as Fig. [TO] but now for 
the total cross-section. 
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